(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     73138,       1987]
NotebookOptionsPosition[     66851,       1775]
NotebookOutlinePosition[     67318,       1793]
CellTagsIndexPosition[     67275,       1790]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"ClebschGordan", "::", "\"\<tri\>\""}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"ClebschGordan", "::", "\"\<phy\>\""}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"SixJSymbol", "::", "\"\<tri\>\""}], "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.456663504453125*^9, 3.45666350465625*^9}, {
  3.5404090219547215`*^9, 3.540409032847345*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
   "oFile", "=", 
    "\"\<C:\\\\Users\\\\work\\\\Data\\\\output\\\\Rb87-D2.txt\>\""}], ";"}], 
  " ", 
  RowBox[{"(*", 
   RowBox[{"output", " ", "file"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
  "pathname", "=", "\"\<C:\\\\Users\\\\work\\\\Data\\\\2012-05-15\\\\\>\""}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"filename", "=", "\"\<Rb-D2-line.dat\>\""}], ";"}], " ", 
  RowBox[{"(*", "datafile", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ExpSpecFile", "=", 
   RowBox[{"pathname", " ", "<>", " ", "filename"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"dataTab", "=", " ", 
    RowBox[{"Import", "[", 
     RowBox[{"ExpSpecFile", ",", "\"\<Table\>\""}], "]"}]}], ";"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "Form", " ", "repetition", " ", "rate", " ", "list", " ", "and", " ", 
    "offset", " ", "frequency", " ", 
    RowBox[{"list", ".", " ", "Extract"}], " ", "from", " ", "actual", " ", 
    "data", " ", "file"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fRepList", "=", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fOffsetList", "=", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n", "=", "1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"While", "[", 
   RowBox[{
    RowBox[{"n", "\[LessEqual]", " ", 
     RowBox[{"Length", "[", "dataTab", "]"}]}], ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"fRepList", "=", 
       RowBox[{"Append", "[", 
        RowBox[{"fRepList", ",", 
         RowBox[{
          RowBox[{"dataTab", "[", 
           RowBox[{"[", "n", "]"}], "]"}], "[", 
          RowBox[{"[", "2", "]"}], "]"}]}], "]"}]}], ";", 
      RowBox[{"(*", 
       RowBox[{"\"\<2\>\"", " ", "=", " ", 
        RowBox[{
        "2", "nd", " ", "column", " ", "of", " ", "the", " ", "datafile", " ",
          "contains", " ", "the", " ", "rep", " ", "rate"}]}], "*)"}], 
      "\[IndentingNewLine]", 
      RowBox[{"fOffsetList", "=", 
       RowBox[{"Append", "[", 
        RowBox[{"fOffsetList", ",", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "2"}], "*", 
           RowBox[{
            RowBox[{"dataTab", "[", 
             RowBox[{"[", "n", "]"}], "]"}], "[", 
            RowBox[{"[", "1", "]"}], "]"}]}], ")"}]}], "]"}]}], ";", " ", 
      RowBox[{"(*", 
       RowBox[{"\"\<1\>\"", " ", "=", " ", 
        RowBox[{
        "1", "st", " ", "column", " ", "of", " ", "the", " ", "data", " ", 
         "file"}]}], "*)"}], "\[IndentingNewLine]", 
      RowBox[{"n", "++"}], ";"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Clear", "[", "n", "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.5407390737883897`*^9, 3.540739088359223*^9}, {
   3.540739992029157*^9, 3.540740041489399*^9}, {3.5407406597673416`*^9, 
   3.5407406605630484`*^9}, {3.540759923421386*^9, 3.54075992699477*^9}, {
   3.540916666913783*^9, 3.5409167107412057`*^9}, {3.5410306226774387`*^9, 
   3.5410306254230385`*^9}, {3.541270497166806*^9, 3.5412704996322856`*^9}, {
   3.541274869409685*^9, 3.541274944053684*^9}, {3.54127535250709*^9, 
   3.5412753550661683`*^9}, {3.5412756166532516`*^9, 3.541275618503356*^9}, 
   3.541296481310985*^9, {3.5413395075363216`*^9, 3.541339527972351*^9}, {
   3.541854954893719*^9, 3.5418550048764396`*^9}, {3.541901835934777*^9, 
   3.5419018359972024`*^9}, {3.5425556289495316`*^9, 3.542555668948188*^9}, {
   3.5425645375795045`*^9, 3.5425645423687353`*^9}, {3.5425708187683387`*^9, 
   3.542570829282873*^9}, {3.5429947403438854`*^9, 3.5429947719235477`*^9}, {
   3.5429967242128005`*^9, 3.54299673992482*^9}, 3.5429987707942295`*^9, {
   3.54860722625*^9, 3.548607242109375*^9}, {3.5486142270302734`*^9, 
   3.5486142323271484`*^9}, {3.5500012398134766`*^9, 3.550001240357422*^9}}],

Cell[CellGroupData[{

Cell["Physical Constants", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell["Speed of Light", "Subsubtitle"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"c", "=", "29979245800."}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"cm", "/", "sec"}], "*)"}]}]], "Input",
 CellChangeTimes->{3.54133953459236*^9}],

Cell[CellGroupData[{

Cell["Planck's Constant", "Subsubtitle"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"hBar", "=", 
    RowBox[{"1.054", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "27"}]]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"ergs", " ", "sec"}], "*)"}]}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Boltzmann Constant", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"kB", "=", 
    RowBox[{"1.380658", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "16"}]]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"erg", "/", "K"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.4566596001875*^9, 3.456659625953125*^9}, {
  3.45665968075*^9, 3.456659686359375*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Atomic Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell["General Parameters", "Subtitle",
 CellChangeTimes->{{3.45666146840625*^9, 3.45666147296875*^9}}],

Cell["Nuclear Spin", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"i", "=", 
   RowBox[{"N", "[", 
    FractionBox["3", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.4730885615*^9, 3.473088561953125*^9}, 3.53632801824225*^9, 
   3.5386062447474585`*^9, {3.541339540382368*^9, 3.541339547522378*^9}}],

Cell[CellGroupData[{

Cell["Transit Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.45666389359375*^9, 3.456663897328125*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[Gamma]T", "=", 
  FractionBox["1", 
   RowBox[{"40.2", " ", 
    SuperscriptBox["10", 
     RowBox[{"-", "9"}]]}]]}]], "Input",
 CellChangeTimes->{{3.456663726421875*^9, 3.4566637388125*^9}, {
   3.45666380203125*^9, 3.456663813921875*^9}, 3.4566922269375*^9, 
   3.5366835360235*^9, 3.538523997111776*^9, {3.539985422708371*^9, 
   3.5399854227239714`*^9}, 3.540776076877905*^9, {3.5407762435859165`*^9, 
   3.540776274178693*^9}, 3.540778024066794*^9, {3.540779180248678*^9, 
   3.540779180857094*^9}, {3.5410999368559904`*^9, 3.54109994055319*^9}, 
   3.541339549542381*^9, {3.5500012637246094`*^9, 3.5500012663085938`*^9}}],

Cell[BoxData["2.487562189054726`*^7"], "Output",
 CellChangeTimes->{3.5500013023671875`*^9, 3.5500013724570312`*^9, 
  3.5500021295878906`*^9, 3.5500022170429688`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"N", "[", "%", "]"}]], "Input",
 CellChangeTimes->{{3.5407780935830674`*^9, 3.540778096063563*^9}}],

Cell[BoxData["2.487562189054726`*^7"], "Output",
 CellChangeTimes->{3.550001302423828*^9, 3.550001372508789*^9, 
  3.5500021296416016`*^9, 3.5500022171103516`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Polarization", "Subsubtitle",
 CellChangeTimes->{{3.456663909203125*^9, 3.456663910671875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"q1", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "1.", ",", "0"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456663912765625*^9, 3.456663917453125*^9}, 
   3.456663957921875*^9, 3.541339552242385*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"q2", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "1.", ",", "0"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{3.4566639629375*^9, 3.5413395529023857`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Parameters", "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659860875*^9}}],

Cell["Ground State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
   3.456659752203125*^9}, 3.45665980753125*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SGnd", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665975521875*^9, 3.456659757859375*^9}, 3.4566598095625*^9, {
   3.540760367051635*^9, 3.5407603674417424`*^9}, {3.5413395550023885`*^9, 
   3.5413395583223934`*^9}}],

Cell[CellGroupData[{

Cell["Ground State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}, {3.540760364617364*^9, 3.5407603656784563`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jGnd", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665975521875*^9, 3.456659757859375*^9}, {3.5413395612023973`*^9, 
  3.541339564522402*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.456662908015625*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinGnd", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jGnd", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxGnd", "=", 
   RowBox[{"jGnd", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Hyperfine A", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
  3.45666040421875*^9}, {3.456660756828125*^9, 3.45666075759375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"aGnd", "=", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"3417341305.452145", ",", "20"}], "]"}]}], " ", ";"}], 
  RowBox[{"(*", 
   RowBox[{"From", " ", "Steck", " ", "review"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
  3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
  3.4566605581875*^9, 3.4566605719375*^9}, {3.456660680015625*^9, 
  3.456660705609375*^9}, {3.45666076115625*^9, 3.4566607615625*^9}, {
  3.47308857578125*^9, 3.473088584546875*^9}, {3.47309640846875*^9, 
  3.473096424140625*^9}, {3.536335670945375*^9, 3.53633567542975*^9}, {
  3.53998552417695*^9, 3.539985574127039*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Hyperfine B", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
   3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
   3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
   3.45666040421875*^9}, 3.4566605775*^9, {3.45666076390625*^9, 
   3.45666076546875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"bGnd", "=", "0."}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
   3.4566605581875*^9, 3.456660602203125*^9}, {3.45666063840625*^9, 
   3.45666063971875*^9}, {3.45666076803125*^9, 3.456660773828125*^9}, 
   3.5413395695324087`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Parameters", "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659887765625*^9}}],

Cell["Intermediate State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}, {3.45665980753125*^9, 3.45665981840625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SInt", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665975521875*^9, 3.456659757859375*^9}, {3.4566598095625*^9, 
  3.45665982178125*^9}, {3.5413395721924124`*^9, 3.5413395762124186`*^9}}],

Cell[CellGroupData[{

Cell["Intermediate State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jInt", "=", 
   RowBox[{"N", "[", 
    FractionBox["3", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.473088588953125*^9, 
   3.47308858946875*^9}, 3.536333582476625*^9, {3.5413396009823647`*^9, 
   3.5413396035823684`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.456662923203125*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinInt", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jInt", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxInt", "=", 
   RowBox[{"jInt", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}, {
  3.456662927125*^9, 3.456662936640625*^9}, {3.540760353210621*^9, 
  3.540760354162483*^9}, {3.5409272522854795`*^9, 3.540927253127901*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Gamma]Int", "=", 
    FractionBox["1", 
     RowBox[{"26.2", " ", 
      SuperscriptBox["10", 
       RowBox[{"-", "9"}]]}]]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Seconds", " ", "From", " ", "Budker"}], ",", "Kimball", ",", " ",
     "DeMille"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.456660078296875*^9, 3.456660092390625*^9}, {3.53668310449225*^9, 
  3.536683106007875*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Center of Gravity Energy", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[Nu]Int", "=", 
   RowBox[{"SetPrecision", "[", 
    RowBox[{"384230484468500", ",", "20"}], "]"}]}], " ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "from", " ", "Steck"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.473088596515625*^9, 3.473088616875*^9}, {
   3.47308883840625*^9, 3.473088843484375*^9}, {3.473096433640625*^9, 
   3.473096439859375*^9}, {3.53633478836725*^9, 3.53633479311725*^9}, 
   3.536607319211*^9, {3.536607447226625*^9, 3.536607498007875*^9}, {
   3.53666646042975*^9, 3.5366664651485*^9}, {3.537893167558856*^9, 
   3.5378931893520937`*^9}, 3.5378943549667425`*^9, 3.5378974743522215`*^9, {
   3.5378985672711487`*^9, 3.53789856825395*^9}, {3.5378998431661897`*^9, 
   3.537899869374236*^9}, 3.5379017571978636`*^9, {3.5386062999866185`*^9, 
   3.538606301137684*^9}, {3.5399856188991175`*^9, 3.539985624343527*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Hyperfine A", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
  3.45666040421875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"aInt", "=", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"84718500", ",", "20"}], "]"}]}], " ", ";"}], 
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "From", " ", "Steck"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
  3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
  3.4566605581875*^9, 3.4566605719375*^9}, {3.473088627078125*^9, 
  3.4730886284375*^9}, {3.473096443796875*^9, 3.473096448328125*^9}, {
  3.5363357313985*^9, 3.536335734726625*^9}, {3.536607079039125*^9, 
  3.53660709955475*^9}, {3.536666472757875*^9, 3.536666476257875*^9}, {
  3.536683146882875*^9, 3.5366831523985*^9}, {3.5372022874831247`*^9, 
  3.5372022919049997`*^9}, {3.5393622605083046`*^9, 3.5393622891187553`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Hyperfine B", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
   3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
   3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
   3.45666040421875*^9}, 3.4566605775*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"bInt", "=", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"12496500", ",", "20"}], "]"}]}], " ", ";"}], 
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "From", " ", "Steck"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
   3.4566605581875*^9, 3.456660602203125*^9}, {3.45666063840625*^9, 
   3.45666063971875*^9}, 3.473088631296875*^9, {3.473096452546875*^9, 
   3.473096456171875*^9}, {3.536335742211*^9, 3.536335745961*^9}, {
   3.536607112711*^9, 3.536607116664125*^9}, 3.53666648311725*^9, {
   3.5366831642735*^9, 3.53668317030475*^9}, {3.5372023036862497`*^9, 
   3.5372023085924997`*^9}, {3.5393622969655685`*^9, 3.5393623013803763`*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Final State Parameters", "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659887765625*^9}, {
  3.456660101109375*^9, 3.4566601018125*^9}}],

Cell["Final State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}, {3.45665980753125*^9, 3.456659829375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SFin", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665975521875*^9, 3.456659757859375*^9}, {3.4566598095625*^9, 
   3.45665983434375*^9}, {3.536335431601625*^9, 3.536335431726625*^9}, 
   3.53634248861725*^9, {3.5413396143423834`*^9, 3.541339617742388*^9}}],

Cell[CellGroupData[{

Cell["Final State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659783546875*^9, 
  3.456659787671875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jFin", "=", 
   RowBox[{"N", "[", 
    FractionBox["5", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.4566597779375*^9, 3.456659780921875*^9}, {3.456661545015625*^9, 
   3.456661545171875*^9}, 3.4569295484212503`*^9, {3.47308863565625*^9, 
   3.473088640640625*^9}, {3.536335439086*^9, 3.536335439226625*^9}, {
   3.541339622032394*^9, 3.5413396261424*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Final State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.45666294725*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinFin", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jFin", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxFin", "=", 
   RowBox[{"jFin", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}, {
  3.456662927125*^9, 3.456662957578125*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Final State Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.456659988796875*^9, 3.45665998953125*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Gamma]Fin", "=", 
    FractionBox["1", 
     RowBox[{"241.", " ", 
      SuperscriptBox["10", 
       RowBox[{"-", "9"}]]}]]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"Seconds", " ", "From", " ", "Heavens"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.4566600681875*^9, 3.45666007421875*^9}, {3.4566616303125*^9, 
  3.456661631765625*^9}, {3.4569295617493753`*^9, 3.4569295620150003`*^9}, {
  3.45951917796875*^9, 3.45951917865625*^9}, {3.536683236757875*^9, 
  3.536683244726625*^9}, {3.538582219701305*^9, 3.5385822378273416`*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Energies, Matrix Elements and Transition Probablities", "Subtitle",
 CellChangeTimes->{{3.456681533796875*^9, 3.456681554234375*^9}}],

Cell["Formula for Calculating Energies with HFS", "Subsubtitle",
 CellChangeTimes->{{3.45666156509375*^9, 3.45666157378125*^9}, {
  3.456662570765625*^9, 3.4566625805625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Kappa]", "[", 
    RowBox[{"f_", ",", "j_"}], "]"}], "=", 
   RowBox[{
    RowBox[{"f", 
     RowBox[{"(", 
      RowBox[{"f", "+", "1"}], ")"}]}], "-", 
    RowBox[{"j", 
     RowBox[{"(", 
      RowBox[{"j", "+", "1"}], ")"}]}], "-", 
    RowBox[{"i", 
     RowBox[{"(", 
      RowBox[{"i", "+", "1"}], ")"}]}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45666169846875*^9, 3.456661713*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Delta]E", "[", 
    RowBox[{"a_", ",", "b_", ",", "f_", ",", "j_"}], "]"}], ":=", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"j", "\[NotEqual]", " ", 
       FractionBox["1", "2"]}], "&&", 
      RowBox[{"i", "\[NotEqual]", 
       FractionBox["1", "2"]}]}], ",", 
     RowBox[{
      RowBox[{
       FractionBox["1", "2"], "a", " ", 
       RowBox[{"\[Kappa]", "[", 
        RowBox[{"f", ",", "j"}], "]"}]}], "+", 
      RowBox[{"b", " ", 
       FractionBox[
        RowBox[{
         RowBox[{
          FractionBox["3", "4"], 
          SuperscriptBox[
           RowBox[{"\[Kappa]", "[", 
            RowBox[{"f", ",", "j"}], "]"}], "2"]}], "+", 
         RowBox[{
          FractionBox["3", "4"], 
          RowBox[{"\[Kappa]", "[", 
           RowBox[{"f", ",", "j"}], "]"}]}], "-", 
         RowBox[{"i", 
          RowBox[{"(", 
           RowBox[{"i", "+", "1"}], ")"}], "j", 
          RowBox[{"(", 
           RowBox[{"j", "+", "1"}], ")"}]}]}], 
        RowBox[{"2", " ", "i", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "i"}], "-", "1"}], ")"}], "j", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "j"}], "-", "1"}], ")"}]}]]}]}], ",", 
     RowBox[{
      FractionBox["1", "2"], "a", " ", 
      RowBox[{"\[Kappa]", "[", 
       RowBox[{"f", ",", "j"}], "]"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456661820734375*^9, 3.45666183171875*^9}, {
   3.456662393875*^9, 3.456662397109375*^9}, {3.456683422046875*^9, 
   3.456683431546875*^9}, 3.460223934290625*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"\[Nu]GndToInt", "[", 
    RowBox[{"f_", ",", "fPrime_"}], "]"}], "=", 
   RowBox[{"\[Nu]Int", "-", 
    RowBox[{"\[Delta]E", "[", 
     RowBox[{"aGnd", ",", "bGnd", ",", "f", ",", "jGnd"}], "]"}], "+", 
    RowBox[{"\[Delta]E", "[", 
     RowBox[{"aInt", ",", "bInt", ",", "fPrime", ",", "jInt"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Nu]GndToFin", "[", 
   RowBox[{"f_", ",", "fDPrime_"}], "]"}], ":=", 
  RowBox[{"2", " ", 
   RowBox[{"Switch", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"f", ",", "fDPrime"}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"2.", ",", "4."}], "}"}], ",", "385284566366300", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"2.", ",", "3."}], "}"}], ",", "385284580777800", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"2.", ",", "2."}], "}"}], ",", "385284592255200", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"2.", ",", "1."}], "}"}], ",", "385284600225200", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"1.", ",", "3."}], "}"}], ",", "385287998122200", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"1.", ",", "2."}], "}"}], ",", "385288009600200", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"1.", ",", "1."}], "}"}], ",", "385288017566900"}], 
    "]"}]}]}]}], "Input",
 CellChangeTimes->{{3.45686233334375*^9, 3.4568623374375*^9}, 
   3.4572661415462503`*^9, {3.47308885534375*^9, 3.473088860125*^9}, {
   3.4730964853125*^9, 3.473096517875*^9}, {3.473097064765625*^9, 
   3.47309707559375*^9}, {3.536683396664125*^9, 3.53668342911725*^9}, {
   3.5385822740464134`*^9, 3.5385822762455387`*^9}, {3.538582433139513*^9, 
   3.538582435955674*^9}, {3.53858286930646*^9, 3.538583146448312*^9}, {
   3.538606352099599*^9, 3.5386064861562667`*^9}, {3.5413396900224895`*^9, 
   3.5413397003025036`*^9}}],

Cell[CellGroupData[{

Cell["Reduced Matrix Elements in J Basis", "Subsubtitle"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"dRedJGndToInt", "=", "1"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"(", "4.51", ")"}], "1.28", " ", 
     SuperscriptBox["10", "6"]}], " ", 
    RowBox[{"(*", 
     RowBox[{"Hz", "/", 
      RowBox[{"(", 
       RowBox[{"V", "/", "cm"}], ")"}]}], "*)"}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"dRedJIntToFin", "=", "1"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     FractionBox["1", 
      RowBox[{"e", " ", "a0"}]], 
     RowBox[{"\[Sqrt]", 
      RowBox[{"(", 
       RowBox[{
        FractionBox[
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", 
             FractionBox["1", "2"]}], "+", "1"}], " ", ")"}], "3", " ", 
          "hBar", " ", 
          SuperscriptBox["c", "3"]}], 
         RowBox[{"4", " ", 
          SuperscriptBox[
           RowBox[{"(", 
            RowBox[{"2", " ", "\[Pi]", " ", "\[Nu]6P12To8S"}], ")"}], "3"]}]],
         "2.04", " ", 
        SuperscriptBox["10", "6"]}], ")"}]}], "1.28", " ", 
     SuperscriptBox["10", "6"]}], " ", 
    RowBox[{"(*", 
     RowBox[{"Hz", "/", 
      RowBox[{"(", 
       RowBox[{"V", "/", "cm"}], ")"}]}], "*)"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.4566056364375*^9, 3.4566056744375*^9}, {
  3.456662686453125*^9, 3.456662724453125*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Reduced Matrix Elements in F Basis", "Subsubtitle",
 CellChangeTimes->{3.4566632656875*^9}],

Cell[BoxData[
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{
       SubscriptBox["dRedGndToInt", 
        RowBox[{"f", ",", "fPrime"}]], "=", 
       RowBox[{"\[Sqrt]", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", "fPrime"}], "+", "1"}], ")"}], 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "f"}], "+", "1"}], ")"}], 
          SuperscriptBox[
           RowBox[{"SixJSymbol", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"jInt", ",", "fPrime", ",", "i"}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{"f", ",", "jGnd", ",", "1"}], "}"}]}], "]"}], "2"], 
          SuperscriptBox["dRedJGndToInt", "2"]}], ")"}]}]}], ";"}], ",", 
     RowBox[{"{", 
      RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], ",",
    "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"f", ",", "fMinGnd", ",", "fMaxGnd"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.456662795671875*^9, 3.456662835546875*^9}, {
  3.456663002453125*^9, 3.456663020171875*^9}, {3.5372063527331247`*^9, 
  3.5372063614518747`*^9}}],

Cell[BoxData[
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      SubscriptBox["dRedIntToFin", 
       RowBox[{"fPrime", ",", "fDPrime"}]], "=", 
      RowBox[{"\[Sqrt]", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", " ", "fDPrime"}], "+", "1"}], ")"}], 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "fPrime"}], "+", "1"}], ")"}], 
         SuperscriptBox[
          RowBox[{"SixJSymbol", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"jFin", ",", "fDPrime", ",", "i"}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{"fPrime", ",", "jInt", ",", "1"}], "}"}]}], "]"}], "2"], 
         SuperscriptBox["dRedJIntToFin", "2"]}], ")"}]}]}], ",", 
     RowBox[{"{", 
      RowBox[{"fDPrime", ",", "fMinFin", ",", "fMaxFin"}], "}"}]}], "]"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.456663051265625*^9, 3.456663093375*^9}, {
   3.456663125046875*^9, 3.456663154296875*^9}, 3.456663192609375*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Matrix Elements in F Basis", "Subsubtitle",
 CellChangeTimes->{{3.456663276328125*^9, 3.456663278140625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Do", "[", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Do", "[", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"Do", "[", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           SubscriptBox["dGndToInt", 
            RowBox[{"f", ",", "fPrime", ",", "m", ",", "mPrime"}]], "=", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"-", "1"}], ")"}], 
             RowBox[{"fPrime", "-", "mPrime"}]], 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"n", "=", "1"}], "3"], 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                "q1", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], 
                RowBox[{"ThreeJSymbol", "[", 
                 RowBox[{
                  RowBox[{"{", 
                   RowBox[{"fPrime", ",", 
                    RowBox[{"-", "mPrime"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"1", ",", 
                    RowBox[{"n", "-", "2"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"f", ",", "m"}], "}"}]}], "]"}]}], ")"}], 
              SubscriptBox["dRedGndToInt", 
               RowBox[{"f", ",", "fPrime"}]]}]}]}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{"m", ",", 
            RowBox[{"-", "f"}], ",", "f"}], "}"}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{"mPrime", ",", 
          RowBox[{"-", "fPrime"}], ",", "fPrime"}], "}"}]}], "]"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], 
    ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"f", ",", "fMinGnd", ",", "fMaxGnd"}], "}"}]}], "]"}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.45666329196875*^9, 3.4566633891875*^9}, {3.456663993*^9,
   3.45666408128125*^9}, {3.45666411909375*^9, 3.456664121515625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Do", "[", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Do", "[", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"Do", "[", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           SubscriptBox["dIntToFin", 
            RowBox[{
            "fPrime", ",", "fDPrime", ",", "mPrime", ",", "mDPrime"}]], "=", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"-", "1"}], ")"}], 
             RowBox[{"fDPrime", "-", "mDPrime"}]], 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"n", "=", "1"}], "3"], 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                "q2", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], 
                RowBox[{"ThreeJSymbol", "[", 
                 RowBox[{
                  RowBox[{"{", 
                   RowBox[{"fDPrime", ",", 
                    RowBox[{"-", "mDPrime"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"1", ",", 
                    RowBox[{"n", "-", "2"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"fPrime", ",", "mPrime"}], "}"}]}], "]"}]}], ")"}], 
              SubscriptBox["dRedIntToFin", 
               RowBox[{"fPrime", ",", "fDPrime"}]]}]}]}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{"mPrime", ",", 
            RowBox[{"-", "fPrime"}], ",", "fPrime"}], "}"}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{"mDPrime", ",", 
          RowBox[{"-", "fDPrime"}], ",", "fDPrime"}], "}"}]}], "]"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"fDPrime", ",", "fMinFin", ",", "fMaxFin"}], "}"}]}], "]"}], 
    ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.45666329196875*^9, 3.456663489203125*^9}, {
  3.456664105703125*^9, 3.45666415221875*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Doppler Distribution Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"amu", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"1.66053886", " ", 
      SuperscriptBox["10", 
       RowBox[{"-", "27"}]]}], ")"}], 
    SuperscriptBox["10", "3"]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mRb", "=", 
   RowBox[{"86.909180520", " ", "amu"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"temp", "=", "323."}], ";"}]}], "Input",
 CellChangeTimes->{{3.456606503203125*^9, 3.4566065051875*^9}, {
   3.456664763890625*^9, 3.456664790515625*^9}, {3.473097352640625*^9, 
   3.473097368078125*^9}, {3.536605457211*^9, 3.53660545749225*^9}, 
   3.5366060920235*^9, {3.536606131101625*^9, 3.536606151414125*^9}, 
   3.53666517880475*^9, {3.540317016135844*^9, 3.540317058781576*^9}, {
   3.540416630451027*^9, 3.5404166312154217`*^9}, {3.5406684546918325`*^9, 
   3.540668459227092*^9}, 3.5406802954125185`*^9, 3.541339635142412*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[Alpha]", "=", 
  RowBox[{"\[Sqrt]", 
   RowBox[{"(", 
    FractionBox[
     RowBox[{"2", " ", "kB", " ", "temp"}], "mRb"], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.45666479665625*^9, 3.456664796984375*^9}, {
  3.536606154070375*^9, 3.536606154336*^9}}],

Cell[BoxData["24860.048218994714`"], "Output",
 CellChangeTimes->{3.5500013041621094`*^9, 3.5500013745126953`*^9, 
  3.5500021320771484`*^9, 3.550002219165039*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"dopplerDist", "[", "v_", "]"}], ":=", 
  RowBox[{"Exp", "[", 
   RowBox[{"-", 
    FractionBox[
     SuperscriptBox["v", "2"], 
     SuperscriptBox["\[Alpha]", "2"]]}], "]"}]}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comb Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell["Offset Frequency", "Subsubtitle"],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"f0", "=", 
    RowBox[{"20.", " ", 
     SuperscriptBox["10", "6"]}]}], ";"}], "*)"}]], "Input",
 CellChangeTimes->{{3.536606229336*^9, 3.536606229476625*^9}, 
   3.536609033945375*^9, {3.53667308924225*^9, 3.536673089507875*^9}, {
   3.536679694632875*^9, 3.5366796956485*^9}, {3.536683666961*^9, 
   3.5366836672735*^9}, 3.5370104677487497`*^9, 3.5370325434049997`*^9, {
   3.539116492685813*^9, 3.5391164930290136`*^9}, 3.539127919839781*^9, 
   3.5400684933016205`*^9, {3.54127557622561*^9, 3.5412755773126698`*^9}}],

Cell[CellGroupData[{

Cell["Rep Rate Range", "Subsubtitle",
 CellChangeTimes->{{3.45666493915625*^9, 3.456664941796875*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"fRepMin", "=", "250000000"}], 
  RowBox[{"(*", 
   RowBox[{"Min", "[", "fRepList", "]"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{"fRepMax", "=", 
  RowBox[{"250000000", "+", "90", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Max", "[", "fRepList", "]"}], ";"}], "*)"}]}]}]}], "Input",
 CellChangeTimes->CompressedData["
1:eJxTTMoPSmViYGCQBWIQzdCQuk7g4EtHBoZ9ELohdBeEL7EPwn9xEMI3PwTh
s56D8OMhdIPTEVkwn/coiG52WHucM/IVkPY7DaYP+LR0gfmL28H0gSbN1WD+
O30IXZe8DkzvSYXQzKEbwOrEI0H0odKbf1WjQPRDKH171how3TwHRL/4/XHz
ymgQ3bAFTFer7V8Npvsh9O+th8+B+fOOgOnfdbP+gGnp2RB66gYIfRpMOz3c
+OsfkC44+hVMf/sZ9h9Ef9n1AUwf4f8hzRzzyvEWo5QMiNapefXYP/aV4xX+
licgWkpilXVI3CtHNqZOexC9YLKAsGrCK8cuTWEREP3kJp+DJpD+c22+I4gu
EeqOOg2kv31dFw2ifXbvOxed+MrRb8UbMG2mtaEiBkgvUt0Nppf9X2KVC6S5
jKXsQHSJlzzDbyD9OCaJCUQDAFIp8ns=
  "]],

Cell[BoxData["250000000"], "Output",
 CellChangeTimes->{3.5500013041845703`*^9, 3.5500013746083984`*^9, 
  3.550002132149414*^9, 3.550002219267578*^9}],

Cell[BoxData["250000090"], "Output",
 CellChangeTimes->{3.5500013041845703`*^9, 3.5500013746083984`*^9, 
  3.550002132149414*^9, 3.5500022192714844`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"fRepAvg", "=", 
   FractionBox[
    RowBox[{"fRepMin", "+", "fRepMax"}], "2"]}], ";"}]], "Input",
 CellChangeTimes->{{3.456665192578125*^9, 3.456665212796875*^9}, {
  3.5407405877168436`*^9, 3.540740594425747*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[Nu]GndToInt", "[", 
  RowBox[{"fMinGnd", ",", "fMinInt"}], "]"}]], "Input",
 CellChangeTimes->{{3.5409273473387165`*^9, 3.540927356371348*^9}}],

Cell[BoxData["3.842344540713818`*^14"], "Output",
 CellChangeTimes->{3.550001304201172*^9, 3.550001374709961*^9, 
  3.5500021322539062`*^9, 3.550002219495117*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comb Modes", "Subsubtitle",
 CellChangeTimes->{{3.4566813895*^9, 3.45668139175*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"n1NomA", "=", 
   RowBox[{"Round", "[", 
    FractionBox[
     RowBox[{"\[Nu]GndToInt", "[", 
      RowBox[{"fMinGnd", ",", "fMinInt"}], "]"}], "fRepAvg"], "]"}]}], 
  RowBox[{"(*", 
   RowBox[{"Round", "[", 
    FractionBox[
     RowBox[{"\[Nu]GndToInt", "[", 
      RowBox[{"1", ",", "1"}], "]"}], "fRepAvg"], "]"}], "*)"}], 
  RowBox[{"(*", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{"n1NomB", "=", 
  RowBox[{"Round", "[", 
   FractionBox[
    RowBox[{"\[Nu]GndToInt", "[", 
     RowBox[{"fMaxGnd", ",", "fMaxInt"}], "]"}], "fRepAvg"], 
   "]"}]}]}], "Input",
 CellChangeTimes->{{3.456665147765625*^9, 3.456665238296875*^9}, {
  3.540759941600395*^9, 3.5407599530539513`*^9}, {3.5407599877579145`*^9, 
  3.540759992454809*^9}, {3.54077629672156*^9, 3.5407763273923397`*^9}}],

Cell[BoxData["1536938"], "Output",
 CellChangeTimes->{3.550001304211914*^9, 3.5500013747978516`*^9, 
  3.550002132786133*^9, 3.5500022195341797`*^9}],

Cell[BoxData["1536912"], "Output",
 CellChangeTimes->{3.550001304211914*^9, 3.5500013747978516`*^9, 
  3.550002132786133*^9, 3.5500022195371094`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"n2NomA", "=", 
   RowBox[{"Round", "[", 
    RowBox[{
     FractionBox["1", "fRepAvg"], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Nu]GndToFin", "[", 
        RowBox[{"fMinGnd", ",", "fMinFin"}], "]"}], "-", 
       RowBox[{"\[Nu]GndToInt", "[", 
        RowBox[{"fMinGnd", ",", "fMinInt"}], "]"}]}], ")"}]}], "]"}]}], 
  RowBox[{"(*", 
   RowBox[{"Round", "[", 
    RowBox[{
     FractionBox["1", "fRepAvg"], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Nu]GndToFin", "[", 
        RowBox[{"1", ",", "3"}], "]"}], "-", 
       RowBox[{"\[Nu]GndToInt", "[", 
        RowBox[{"2", ",", "1"}], "]"}]}], ")"}]}], "]"}], "*)"}], 
  RowBox[{"(*", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{"n2NomB", "=", 
  RowBox[{"Round", "[", 
   RowBox[{
    FractionBox["1", "fRepAvg"], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"\[Nu]GndToFin", "[", 
       RowBox[{"fMaxGnd", ",", "fMaxFin"}], "]"}], "-", 
      RowBox[{"\[Nu]GndToInt", "[", 
       RowBox[{"fMaxGnd", ",", "fMaxInt"}], "]"}]}], ")"}]}], "]"}], 
  RowBox[{"(*", 
   RowBox[{"Round", "[", 
    RowBox[{
     FractionBox["1", "fRepAvg"], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Nu]GndToFin", "[", 
        RowBox[{"2", ",", "4"}], "]"}], "-", 
       RowBox[{"\[Nu]GndToInt", "[", 
        RowBox[{"1", ",", "1"}], "]"}]}], ")"}]}], "]"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.456665147765625*^9, 3.4566653128125*^9}, {
  3.5407600090265756`*^9, 3.5407600152370872`*^9}, {3.5407600962390084`*^9, 
  3.5407601177729425`*^9}, {3.5407602283294077`*^9, 3.5407602651555557`*^9}, {
  3.540776333772985*^9, 3.5407764024780273`*^9}}],

Cell[BoxData["1545366"], "Output",
 CellChangeTimes->{3.550001304225586*^9, 3.5500013748984375`*^9, 
  3.5500021329257812`*^9, 3.550002219571289*^9}],

Cell[BoxData["1545364"], "Output",
 CellChangeTimes->{3.550001304225586*^9, 3.5500013748984375`*^9, 
  3.5500021329257812`*^9, 3.5500022195751953`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"n1List", "=", 
  RowBox[{"Union", "[", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n1NomA", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n1NomB", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "]"}]}], "]"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.456665383359375*^9, 3.45666546196875*^9}, {
   3.456665618125*^9, 3.456665620609375*^9}, {3.5372066701081247`*^9, 
   3.5372066730768747`*^9}, {3.5377187815151253`*^9, 
   3.5377187956088753`*^9}, {3.5377946040932503`*^9, 
   3.5377946077338753`*^9}, {3.5377949385776253`*^9, 
   3.5377949450151253`*^9}, {3.537893207650926*^9, 3.537893217432143*^9}, {
   3.5385234618671618`*^9, 3.53852346744048*^9}, {3.538583211301021*^9, 
   3.5385832213905983`*^9}, {3.54031706878374*^9, 3.5403170739486637`*^9}, {
   3.540409056767713*^9, 3.5404090620360146`*^9}, {3.5404810791778603`*^9, 
   3.5404810845130253`*^9}, 3.5407000000628*^9, 3.540740600775801*^9, {
   3.5407601382145753`*^9, 3.54076014211565*^9}, {3.5407763928212557`*^9, 
   3.5407763961909857`*^9}, {3.540776760262188*^9, 3.5407767670640497`*^9}, {
   3.5407771653785686`*^9, 3.540777170464364*^9}, {3.540917036115741*^9, 
   3.540917040328416*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1536909", ",", "1536910", ",", "1536911", ",", "1536912", ",", "1536913", 
   ",", "1536914", ",", "1536915", ",", "1536935", ",", "1536936", ",", 
   "1536937", ",", "1536938", ",", "1536939", ",", "1536940", ",", 
   "1536941"}], "}"}]], "Output",
 CellChangeTimes->{3.5500013042373047`*^9, 3.5500013753828125`*^9, 
  3.5500021330527344`*^9, 3.5500022196015625`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"n2List", "=", 
  RowBox[{"Union", "[", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n2NomA", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n2NomB", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "]"}]}], "]"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.456665383359375*^9, 3.45666546196875*^9}, {
   3.456665585125*^9, 3.45666558903125*^9}, {3.456665624328125*^9, 
   3.456665626734375*^9}, {3.5372066753268747`*^9, 3.5372066804518747`*^9}, {
   3.5377187866401253`*^9, 3.5377187978745003`*^9}, {3.5377946088901253`*^9, 
   3.5377946116245003`*^9}, {3.5377949469370003`*^9, 
   3.5377949552338753`*^9}, {3.5378932198813477`*^9, 3.537893225029357*^9}, {
   3.538523470201638*^9, 3.5385234737538414`*^9}, {3.538583224257762*^9, 
   3.5385832272179317`*^9}, {3.5403170767417803`*^9, 3.540317081079692*^9}, {
   3.5404090654172077`*^9, 3.540409069285429*^9}, {3.5404810860886154`*^9, 
   3.5404810923285756`*^9}, 3.540699997145656*^9, 3.540740602460828*^9, {
   3.540760133049552*^9, 3.540760136685354*^9}, {3.5407767689985237`*^9, 
   3.5407767765648146`*^9}, {3.5407771721336284`*^9, 
   3.5407771762989883`*^9}, {3.540917041857461*^9, 3.5409170471155033`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1545360", ",", "1545361", ",", "1545362", ",", "1545363", ",", "1545364", 
   ",", "1545365", ",", "1545366", ",", "1545367", ",", "1545368", ",", 
   "1545369", ",", "1545370"}], "}"}]], "Output",
 CellChangeTimes->{3.5500013042470703`*^9, 3.550001375873047*^9, 
  3.5500021332753906`*^9, 3.550002219631836*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comb Equation", "Subsubtitle",
 CellChangeTimes->{{3.456681410328125*^9, 3.456681413015625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Omega]L", "[", 
    RowBox[{"n_", ",", "fR_", ",", "f0_"}], "]"}], ":=", 
   RowBox[{"2", " ", "\[Pi]", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"n", " ", "fR"}], "+", "f0"}], ")"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456678515171875*^9, 3.45667856759375*^9}, {
  3.4566795630625*^9, 3.45667956609375*^9}, {3.541275345195319*^9, 
  3.541275369178735*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Integrate 2-photon Lineshape over Doppler Distribution", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"transProb", "[", 
     RowBox[{"\[Omega]1_", ",", "\[Omega]2_", ",", "v_"}], "]"}], "=", 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"fDPrime", "=", "fMinFin"}], "fMaxFin"], 
     RowBox[{
      UnderoverscriptBox["\[Sum]", 
       RowBox[{"f", "=", "fMinGnd"}], "fMaxGnd"], 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         FractionBox[
          SuperscriptBox[
           RowBox[{"(", 
            RowBox[{"\[Gamma]Fin", "+", "\[Gamma]T"}], ")"}], "2"], "4"], "/", 
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{"2", " ", "\[Pi]", " ", 
               RowBox[{"\[Nu]GndToFin", "[", 
                RowBox[{"f", ",", "fDPrime"}], "]"}]}], "-", 
              RowBox[{
               RowBox[{"(", "\[Omega]1", ")"}], 
               RowBox[{"(", 
                RowBox[{"1", "+", 
                 FractionBox["v", "c"]}], ")"}]}], "-", 
              RowBox[{
               RowBox[{"(", "\[Omega]2", ")"}], 
               RowBox[{"(", 
                RowBox[{"1", "-", 
                 FractionBox["v", "c"]}], ")"}]}]}], ")"}], "2"], "+", 
           FractionBox[
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"\[Gamma]Fin", "+", "\[Gamma]T"}], ")"}], "2"], "4"]}], 
          ")"}]}], "\[IndentingNewLine]", 
        FractionBox["1", 
         RowBox[{
          RowBox[{"2", "f"}], "+", "1"}]], 
        RowBox[{
         UnderoverscriptBox["\[Sum]", 
          RowBox[{"m", "=", 
           RowBox[{"-", "f"}]}], "f"], 
         RowBox[{
          UnderoverscriptBox["\[Sum]", 
           RowBox[{"mDPrime", "=", 
            RowBox[{"-", "fDPrime"}]}], "fDPrime"], 
          SuperscriptBox[
           RowBox[{"Abs", "[", 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"fPrime", "=", "fMinInt"}], "fMaxInt"], 
             RowBox[{
              UnderoverscriptBox["\[Sum]", 
               RowBox[{"mPrime", "=", 
                RowBox[{"-", "fPrime"}]}], "fPrime"], "\[IndentingNewLine]", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 FractionBox[
                  RowBox[{"\[Gamma]Int", "+", "\[Gamma]T"}], "2"], "/", 
                 RowBox[{"(", 
                  RowBox[{
                   FractionBox[
                    RowBox[{"\[ImaginaryI]", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Gamma]Int", "+", "\[Gamma]T"}], ")"}]}], "2"], 
                   "+", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"(", "\[Omega]1", ")"}], 
                    RowBox[{"(", 
                    RowBox[{"1", "+", 
                    FractionBox["v", "c"]}], ")"}]}], "-", " ", 
                    RowBox[{"2", " ", "\[Pi]", " ", 
                    RowBox[{"\[Nu]GndToInt", "[", 
                    RowBox[{"f", ",", "fPrime"}], "]"}]}]}], ")"}]}], ")"}]}],
                 ")"}], 
               SubscriptBox["dGndToInt", 
                RowBox[{"f", ",", "fPrime", ",", "m", ",", "mPrime"}]], 
               SubscriptBox["dIntToFin", 
                RowBox[{
                "fPrime", ",", "fDPrime", ",", "mPrime", ",", 
                 "mDPrime"}]]}]}]}], "]"}], "2"]}]}]}], ")"}]}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"fnToIntegrate", "[", 
     RowBox[{"fR_", ",", "fOffset_", ",", "v_"}], "]"}], "=", 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"n1", "=", "1"}], 
      RowBox[{"Length", "[", "n1List", "]"}]], 
     RowBox[{
      UnderoverscriptBox["\[Sum]", 
       RowBox[{"n2", "=", "1"}], 
       RowBox[{"Length", "[", "n2List", "]"}]], 
      RowBox[{
       RowBox[{"dopplerDist", "[", "v", "]"}], 
       RowBox[{"transProb", "[", 
        RowBox[{
         RowBox[{"\[Omega]L", "[", 
          RowBox[{
           RowBox[{
           "n1List", "\[LeftDoubleBracket]", "n1", "\[RightDoubleBracket]"}], 
           ",", "fR", ",", "fOffset"}], "]"}], ",", 
         RowBox[{"\[Omega]L", "[", 
          RowBox[{
           RowBox[{
           "n2List", "\[LeftDoubleBracket]", "n2", "\[RightDoubleBracket]"}], 
           ",", "fR", ",", "fOffset"}], "]"}], ",", "v"}], "]"}]}]}]}]}], 
   ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vIntSpec", "[", 
    RowBox[{
     RowBox[{"fRep_", "?", "NumericQ"}], ",", 
     RowBox[{"fO_", "?", "NumericQ"}]}], "]"}], ":=", 
   RowBox[{"NIntegrate", "[", 
    RowBox[{
     RowBox[{"fnToIntegrate", "[", 
      RowBox[{"fRep", ",", "fO", ",", "v"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"v", ",", 
       RowBox[{"-", "60000"}], ",", "60000"}], "}"}], ",", 
     RowBox[{"MinRecursion", "\[Rule]", "3"}], ",", 
     RowBox[{"MaxRecursion", "\[Rule]", "10"}], ",", 
     RowBox[{"AccuracyGoal", "\[Rule]", "4"}], ",", 
     RowBox[{"PrecisionGoal", "\[Rule]", "2"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.538523548469115*^9, 3.5385235524603434`*^9}, {
  3.5385832601238136`*^9, 3.538583261230877*^9}, {3.54077684144771*^9, 
  3.5407768433509836`*^9}, {3.5407781956265917`*^9, 3.5407781988247147`*^9}, {
  3.540778342428238*^9, 3.540778354378298*^9}, {3.54090952236373*^9, 
  3.5409095316782627`*^9}, {3.5409167735568705`*^9, 3.5409167781128006`*^9}, {
  3.5409170292350388`*^9, 3.5409170310449286`*^9}, {3.5409218751219053`*^9, 
  3.5409218761359835`*^9}, {3.54127538747349*^9, 3.5412754413861923`*^9}, {
  3.5412754775292873`*^9, 3.5412754971543236`*^9}, {3.541854908873424*^9, 
  3.541854930573163*^9}, {3.541901821732953*^9, 3.5419018296297913`*^9}, {
  3.5486142020615234`*^9, 3.5486142102021484`*^9}, {3.5486152438740234`*^9, 
  3.5486152532802734`*^9}, {3.5486153189365234`*^9, 3.5486153292021484`*^9}, {
  3.5486156190771484`*^9, 3.5486156448271484`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Calculation", "Subtitle",
 CellChangeTimes->{{3.456681446140625*^9, 3.45668144753125*^9}, {
  3.5379007941050625`*^9, 3.537900798020669*^9}},
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"\[Delta]f", "=", "0.25"}], ";", " ", 
   RowBox[{"(*", 
    RowBox[{"step", " ", "size", " ", "in", " ", "rep", " ", "rate"}], "*)"}],
    "\[IndentingNewLine]", 
   RowBox[{"fr", "=", "250000000"}], ";", "\[IndentingNewLine]", 
   RowBox[{"fRepList", "=", 
    RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"specList", "=", 
    RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"Form", " ", "the", " ", "fRepList"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{"While", " ", "[", 
    RowBox[{
     RowBox[{"fr", "<", "250000090"}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"fRepList", "=", 
        RowBox[{"Append", "[", 
         RowBox[{"fRepList", ",", "fr"}], "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"fr", "+=", "\[Delta]f"}], ";"}], "}"}]}], "]"}]}], 
  "*)"}]], "Input",
 CellChangeTimes->{
  3.456746939*^9, {3.5372003913424997`*^9, 3.5372003937643747`*^9}, {
   3.5378885919057503`*^9, 3.5378885935932503`*^9}, {3.537901719773398*^9, 
   3.53790172663741*^9}, 3.5379017983975363`*^9, 3.5385235615778646`*^9, 
   3.538524008413422*^9, {3.540084146036523*^9, 3.5400841470661955`*^9}, 
   3.5404312983500743`*^9, 3.5404424331582637`*^9, {3.5406998704760923`*^9, 
   3.540699870788086*^9}, {3.5407000208571997`*^9, 3.5407000214031897`*^9}, {
   3.5407003545674286`*^9, 3.540700373396387*^9}, {3.5407005292071896`*^9, 
   3.54070054245142*^9}, {3.540701011319009*^9, 3.5407010167789383`*^9}, {
   3.540701078148552*^9, 3.54070107905334*^9}, {3.540702894168069*^9, 
   3.54070290344995*^9}, 3.540704212332471*^9, 3.54070429143855*^9, {
   3.5407339187505383`*^9, 3.5407339258319435`*^9}, {3.5407346581998324`*^9, 
   3.5407346626740885`*^9}, {3.5407390969557147`*^9, 3.540739098918827*^9}, 
   3.540740429510459*^9, {3.5410306331138387`*^9, 3.5410306802726383`*^9}, {
   3.5410308040274386`*^9, 3.541030814635439*^9}, {3.5412704895519075`*^9, 
   3.5412704916584883`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"rFrep", "=", "0.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"tt", "=", 
    RowBox[{"AbsoluteTime", "[", "]"}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
   "distribute", " ", "NIntegrate", " ", "among", " ", "CPU", " ", "cores"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"specList", "=", 
   RowBox[{"{", "}"}]}], ";", 
  RowBox[{
  "SetSharedVariable", "[", "specList", "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ParallelDo", "[", 
   RowBox[{
    RowBox[{"AppendTo", "[", 
     RowBox[{"specList", ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
         "fRepList", "\[LeftDoubleBracket]", "kk", "\[RightDoubleBracket]"}], 
         ")"}], ",", 
        FractionBox[
         RowBox[{"Chop", "[", 
          RowBox[{"vIntSpec", "[", 
           RowBox[{
            RowBox[{
            "fRepList", "\[LeftDoubleBracket]", "kk", 
             "\[RightDoubleBracket]"}], ",", 
            RowBox[{
            "fOffsetList", "\[LeftDoubleBracket]", "kk", 
             "\[RightDoubleBracket]"}]}], "]"}], "]"}], "12"]}], "}"}]}], 
     "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"kk", ",", 
      RowBox[{"Length", "[", "fRepList", "]"}]}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ttt", "=", 
   RowBox[{
    RowBox[{"AbsoluteTime", "[", "]"}], "-", "tt"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{"\"\<Execution time = \>\"", ",", "ttt"}], "]"}]}], "Input",
 CellChangeTimes->{{3.540700569251876*^9, 3.540700631292281*^9}, {
   3.540704216154398*^9, 3.540704252064907*^9}, {3.5407042969296446`*^9, 
   3.5407042987392097`*^9}, 3.5407055725228095`*^9, {3.54073392927214*^9, 
   3.540734037008302*^9}, {3.541275778337726*^9, 3.541275794685014*^9}, {
   3.541276893218398*^9, 3.541276895647437*^9}, {3.541296195763912*^9, 
   3.541296442665784*^9}, {3.5412965326870537`*^9, 3.5412965488659134`*^9}, {
   3.541298036226572*^9, 3.5412980663070354`*^9}, {3.5412999438710012`*^9, 
   3.5413000247309847`*^9}, 3.5413009553674545`*^9, {3.541302060074316*^9, 
   3.54130207288249*^9}, {3.541302810233549*^9, 3.5413028515127363`*^9}, {
   3.541304182880085*^9, 3.54130419865219*^9}, {3.542556980744997*^9, 
   3.542556983272213*^9}, {3.5425636106215625`*^9, 3.542563611370367*^9}, {
   3.5500013545*^9, 3.550001365375*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Execution time = \"\>", "\[InvisibleSpace]", 
   "2600.2724609`10.866563849953662"}],
  SequenceForm["Execution time = ", 2600.2724609`10.866563849953662],
  Editable->False]], "Print",
 CellChangeTimes->{3.5500048204628906`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"totalSpec", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{
         "fRepList", "\[LeftDoubleBracket]", "t", "\[RightDoubleBracket]"}], 
         ",", 
         RowBox[{"Chop", "[", 
          RowBox[{
          "specList", "\[LeftDoubleBracket]", "t", "\[RightDoubleBracket]"}], 
          "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "1", ",", 
         RowBox[{"Length", "[", "specList", "]"}]}], "}"}]}], "]"}]}], ";", 
    "\n", 
    RowBox[{"fRepList", "=", 
     RowBox[{"fRepList", "-", 
      RowBox[{"Min", "[", "fRepList", "]"}]}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"plotSpec", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{
         "fRepList", "\[LeftDoubleBracket]", "t", "\[RightDoubleBracket]"}], 
         ",", 
         RowBox[{"Chop", "[", 
          RowBox[{
          "specList", "\[LeftDoubleBracket]", "t", "\[RightDoubleBracket]"}], 
          "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "1", ",", 
         RowBox[{"Length", "[", "specList", "]"}]}], "}"}]}], "]"}]}], ";"}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{"ListPlot", "[", 
   RowBox[{
    RowBox[{"Sort", "[", "specList", "]"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
    RowBox[{"Joined", "\[Rule]", "True"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5409271816468678`*^9, 3.5409272231595325`*^9}, 
   3.540959220519678*^9, {3.541299917754927*^9, 3.5412999284104104`*^9}, {
   3.5413019390284843`*^9, 3.5413019427882533`*^9}, {3.5413033228068624`*^9, 
   3.5413033265042048`*^9}, {3.541303499156045*^9, 3.5413034998112698`*^9}}],

Cell[BoxData[
 GraphicsBox[
  {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwl2HlYTVsDBvAtpDKUXBSZCWmSeagUhVBoMCRChZCM6RYdV2mk6TRSKSpp
LikiCblSUjJFMpSUVCrXrK/3/f7q+T1va9xrr7XOHrNlz2pbCUEQtvcQBPx1
KujqsilL1Ru62l7KICBTVxBEgm23hR3hqusHZujeGS30gI/1aPkrSEjTnVfr
RqtMPz38ZnxKt7tozR6v5RPrk5FLwBVeK45rXk3urk/oCSuuarVaNyaluz4R
vUNmjZN4dirq7wXnnJw2PEQuHe3TegqHIsobMlC+N2w90r3ujpCNnD5z00Sq
54zLyCVZ38MSiyOFV5HTxyJrxNEJhegfHb3h+uaRKneQ9+H4zpUrSybf1+0e
P932wPlr3d8VyKXgAT1W7NmV8xj9k4ZL25/4P49/hpzWbj9e5uPwHLkMLB2V
u7hG/yly2n5HW271u0foX194royvU+7acoyfnvBK5Ytf9b/oHy1rPL2vQnUR
8n7wHwnxj09e11F/f9jRRU36g08u5pdev6ktaOUGzIcwAG7RkXBomp6hmxXl
Rq84Hby+9EQKclk4ckdV2Zh9SbreeiI6Y9T34PFvEpDLcbwmlTdqpRPQPq3x
9dvX/ePi0T/6YeLGd2lp8RjfQNggJHju9EKWl4dnSuivGZZzATl95pb68buT
+fwHwQtcNuS7x6ai/7Rhx9LABVPxvIW/4Bnrtlsty0hDefrgr2OB7bZYL8Jg
OLBmYeE5w0TkdIz/+5X2w9EfYQjrL99fUt8vDjltvTckKXBaDPpPxx3t8m4N
P418KLxZRtDaOigCOR2ae2NGin4Yxq/A+v7rnb2xJQTtK8IhN2Ltv7SKkdMq
E93HZ2nBwjD4gewipd4JwcjpuDi18zOnBKN/w+Gfy9+KVDuDMH76gbfJVu+S
IJRXghd/ttJujWdOy9cpLKnJDUL5EfDftTlpXf7MaUFWXGeaCnfRPz1mGNdW
4v9FI/k+OlzyLpFh+6Ngd9VVus4r2D+6QNd+sJFrMNY/vdt6w6ehv4NR/2h4
rfvE+iNDxchpR5mKvgoNYtQ/Br7qJMSmjOX4x8JKSs2RzpNCUT+tnqrortEY
ivbHwV8i5BLv3OD80huSwhqDMpiPh/tP7y2VdCcCOa1Q/UDBojgC+QR4u7TR
bfuvp9E/eoaUslzWoyi0rwyXdTWoTAyPQf9o0VSzooWlZzE/9NL2cg+X9tju
8bhNhG/vSy0qMDiH8pOYv/T2dvU8j/L00XYzL8kjCShP5zU2P+1pewHtT4bV
8lMiJAyx33XRj69GTbwyn+tXBT67X7PgvGQm6p8Cd3Qsavc0ycb46KrmravT
O3MwPlXuL/YrvuX5XkF5+o6KnXHd+ALkaty/PL+pH9pchPbp3vctZ21bW4z6
1dnfCqf8XsGlGB8tm266d+zRSpTXgEe2Fb/Y4sT9iw6w8t7jt+Q1ck2+v9l+
sfKr65HT7RL/fBjr3oR8KtfjqxFpeWktyGkJ1Q8tE2taMX66rqUu8OoW5lpc
D5uzRC6Dm5DT+q/KpYY9q0M+Da4coew9cQDbn8731WjLd/++zzA+ukP55uTz
qez/DDhCveDicMcy5LTF3ITbFp13UT994lZw2/iRtzG/M3nevFD9oGR6E+Vn
cb/7/SLG5hvmU0RbSXsmaepfQ3naWdVUzs6Z58dsOPOEuozY9gpyurSySOp3
ex7anwOfVHJp92rLQ/1z4foJ/xXIr81D+/SX0LZH4+bnI5/H/XeEz6Zx3ft5
d/30tlGf5/ddV9i9HwvzYbPP+X1Myrn/03MyTAaPUS5GeW1YbtE9i31S95HT
fpd0r1dJ4rwSdPi8mqaYV3x8jPrp1eaaTc3BL7E+dOH/lk8yb9B7i/K03sbZ
t8ateI98AddHmNrR11MbUZ5ectbgz7onjRg/rfOqXcfBpQG5HtvTH6PmehrP
U9CHixxFW+u98DxF9L3L6vOmPH2OfCE83rttufzSKrRPDxuZ+6RfxEPki+BX
ria16atKUZ62u62Y+GH9PbRPGx+W00isuovyBlzvF3t9zA2DBUP4lNZ/P5UC
/sX6px09ylKUbpQgX8z1YDqvqSrhAcrTB/Lf1cQ5ViJfwvlZlpKrLPcUOf07
fo2G+94a5Eu5ny9p62Ol+w45ndGgLcjafUD/jbieXH92PlD4hPVBT92Q9VDm
cSvyZdwvSiP2eZrgfXGjNbfFi8b2aUa+nO+T67Jju+U5v/TXAbqp2nfeoP0V
/9//vnTNk3yB9mkdw1mvExwxn4Ixn2dBVPcGV47y9FiNyttfU/5F/SZwD0un
KI2nt1CeltO751prVYh8JfxslpWei2s++kcr3lU7US7P9b0K9niS+GXC92zk
9IUR2y6G/8Z9UVjN/fde9ti+UrivieggqwnLx93g/cOU/ZH43ivoZxLap6/b
7Ty5PjMRz5eucp0vbXYjHrkZXKiSG76jL/ZjwRzW9nQONBTHIacf9z/hn27P
/Z229TodmzKU+78F/NzRLvblhmjk9OjjEaKD4jPI13D8Ve6XGroiMX9ruV7O
2KicScZ5I6LnBMzbedg2HPk6vu9D3KZUG4Zh/PT63R+6LOaG4nykX/1V2yJe
GYLy69le1O8Jnlo4LwVLPr8Ng4qjQnAei2jj+UHR/ZR4Pm+A38+X2L3dBfcB
EX0tyj3JpiEQuRXXy/on8a3WgcjpnRb62b71ARgfXaC1e7TCgQDUvxF+s+rA
3ExJWNgE7+of4dPD0h/9p48pn3o7JesUcmv4z/HDZ6/JnUJOf9Yy0H7kcBLj
3wxPMl8z/FarH+qnpSumDvFx9sP+tQUebODxPqOvH8rTSsU5Fkkxvqh/K9f/
t6PNygtgET1hwTa5MC8f1G8DRw1yjwnc5IP3l96r4hydauKD8rbc35t7O+XN
gUX08x67oiT6M7eDf3nmDx/a4I326SMf7evU7nqj/m08P0ZvboxM98b9lz60
M6T6UIQ3nh99rMh7koGY5bfDknMMPGwCvVH/DpYf2H9CUxhzujO6c61VKHN7
2D5/0falR2AR/b6Xqk5sOPOd8IC2vypEDczpb2E7RwVVMN/F/eRbUfqOBZgP
ET14j37kYQ3Oz27YoevWvFd5Pmifbrp2xfDGD7iLlvhT+XbgSl+Ud4DPVt2d
eeWCL8ZHd77Uq73Uxw/5HjjeYuq0Zjs8T8ER3uczOrLpCZ8ffVJ2Snm62UnU
T/sZTq/VfHoS/d/L/mSV1yxcwvWzD57dmhQddgnrRURXH3d6dP2yP/L9XL++
HYWL5mO9utEJW/Pe/Urm+jwA112v7lgnH4icNvL09OsUc/0fhLuuSD05IcX3
g7bpqL9ZnMT7Lz04cI3JFqNgPN9DcMPCezWW33DfFZy4/t+ftr6TIUb99LfS
/hnZ+nhfhcNw+M+ynh598D670cKD2E0j23CfFZx5PvQ8knRODfuBiB5c+6c+
fBHvp3/DOeFP2q66RqI8Hb6warpzPu6rXfQI8dx2hTm4r4pcOB+Pt9b9k839
iT5V8sVh3vazmF9XOMzDwMG6ERaOwNFrqlZZdcYipxX3b7LZO+sc8qPwVSNV
84WLuX/SP55lFk7r4P2VLq12UzxacQH9c4N3940MetF0EeVFcIiidPTk9amo
n1Y7cPB1pFIGytOrwoqlvW/x+8AxrlfLuNaIXJwHXXSJ5HeTC9L4PS76Bx7z
WTIkz5Tnx3HOV1FEhbHONd2n50R0xzbdMZHz+fucvve48LvZoyKUd+f81Q8v
UppSjJw+s0xTR04V9wE3D/jB59je1mkPUP8JrvfjJnJL5Hm/oMNdbKZ9DeX9
w5PzUXdQembqa+R0UYhJoaT6e9RPH5RsTzNZ/Qnte8GzDPrpVgZ1oLw3z4ul
fwIiYr6hfXq5gWhA3ItfyH3gGaoWs7Rj/6A8vTfF61HLot9YP75wf2m1+8km
35HT8YW/dJpXdGL/9IOvyPTp1G1qQf20mmVH+nUl3AeEk1wvp5xrlq17g/K0
QdTYkU0h1eg/7SderDqitQr5KZ6f/WpmrX7P+5U/98+hiZlRb+6jflprR0Ce
nOJd5AHwQ8v3iY6lt1CeTpj84cdODd6fA3l+r7157fhB3p9px4wLx5ONeL8N
giunJhiNK+b9mA7376lZd4m/b4JhnTZl2Z7bYBHtPkT13YAxl5CL4QUZUpET
5vH3Eb38bsnsIq8szF8IbJ6zeUGUchbqpw+42YnNh2ehfCh883zlYUVtWESb
DfSY+KyYeRgcdH2/VFs+vx/RfSTKVolkc1B/OCy1YrKS0uZctE8fbh60cKj7
VeQRsJ2JvWGj7A2Up11/hWp+fIr1KkTy/B+o4VuSgN8jbrT7uY9qBXse4PnQ
S2pWrht1/hHqP818t2nPzHT83umiPZbfrzkfg/um6AzcT0dsq2z0GusjCm43
+6D1u6QW7dM/Xti1vNGvQfvRzC0N3C8YPcP7TWvLJ67p3R/tCTGwWG9vwubJ
ZaifLoiemZZ/gO8X/XJ01f5xV/G8RWfhT6qH/ri/vYbzgzZOCIt1U+Hvm1j+
f6ZKo0EFn18c9y9L/2kvt3P+6bu+qg5lyfx+eA7e8vxQ3MI3/H1MyzjUX45W
x/egLnrlI//Zv9RxXxSdhy3GFavdupOsmxXVRdfleEQad/H7Yjxco26dHHCH
34cS4Mw8P11lPdZPL9IO/Soy5ffDRFjhmvW+1yOwnkR0eOLtOJfiSxgfrfks
1O2gmOv3ArxdIyr8iTJ/fyexvfKOmuihfD9orUvN2e7X+f2OTug8sFH6Zznm
/yLXe43N8x/SVcjpgXKmqdNDnqL+ZO4X5r7VI3Zxf0qBvcxs38ec5+9fet2W
0M9KEfg9J6TCubc3Hpz4kt8n6cJNzj2DL99Hnsb9Qr7fEmUFft+k7XaGKm7R
4v0+nfPv9/l1P29+X6AtJbIL1Pbw+WXAv0JzUjs2cX+n7yVXDUtoxHkgZMKr
591bqb0f54WI7p3r/3FhH3zvE7LY/9NzIwYZ8/shrWKuVdL2D84jIRsef/nx
vWhrfj+k64ZkmdfIJuj+D9L+ScI=
    "]]},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  PlotRangeClipping->True]], "Output",
 CellChangeTimes->{3.5500048206191406`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Export", "[", 
  RowBox[{"oFile", ",", 
   RowBox[{"Sort", "[", "specList", "]"}], ",", "\"\<Table\>\""}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.541303835300945*^9, 3.5413038495754023`*^9}}],

Cell[BoxData["\<\"C:\\\\Users\\\\work\\\\Data\\\\output\\\\Rb87-D2.txt\"\>"], \
"Output",
 CellChangeTimes->{3.5500048207753906`*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1272, 684},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
CellContext->Notebook,
Magnification:>FEPrivate`If[
  FEPrivate`Equal[FEPrivate`$VersionNumber, 6.], 1.25, 1.25 Inherited],
FrontEndVersion->"8.0 for Microsoft Windows (32-bit) (November 7, 2010)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 486, 13, 80, "Input"],
Cell[1034, 35, 3933, 94, 322, "Input"],
Cell[CellGroupData[{
Cell[4992, 133, 76, 1, 62, "Subtitle"],
Cell[5071, 136, 37, 0, 35, "Subsubtitle"],
Cell[5111, 138, 190, 6, 35, "Input"],
Cell[CellGroupData[{
Cell[5326, 148, 40, 0, 35, "Subsubtitle"],
Cell[5369, 150, 223, 8, 37, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5629, 163, 107, 1, 35, "Subsubtitle"],
Cell[5739, 166, 330, 10, 37, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[6118, 182, 75, 1, 62, "Subtitle"],
Cell[6196, 185, 102, 1, 62, "Subtitle"],
Cell[6301, 188, 149, 2, 35, "Subsubtitle"],
Cell[6453, 192, 323, 7, 53, "Input"],
Cell[CellGroupData[{
Cell[6801, 203, 105, 1, 35, "Subsubtitle"],
Cell[CellGroupData[{
Cell[6931, 208, 653, 12, 53, "Input"],
Cell[7587, 222, 167, 2, 35, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7791, 229, 122, 2, 35, "Input"],
Cell[7916, 233, 163, 2, 35, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[8128, 241, 101, 1, 35, "Subsubtitle"],
Cell[8232, 244, 242, 6, 35, "Input"],
Cell[8477, 252, 192, 5, 35, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[8718, 263, 106, 1, 62, "Subtitle"],
Cell[8827, 266, 222, 3, 35, "Subsubtitle"],
Cell[9052, 271, 361, 8, 53, "Input"],
Cell[CellGroupData[{
Cell[9438, 283, 247, 3, 35, "Subsubtitle"],
Cell[9688, 288, 285, 7, 53, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10010, 300, 109, 1, 35, "Subsubtitle"],
Cell[10122, 303, 301, 8, 58, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10460, 316, 389, 5, 35, "Subsubtitle"],
Cell[10852, 323, 920, 17, 35, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11809, 345, 413, 6, 35, "Subsubtitle"],
Cell[12225, 353, 589, 10, 35, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[12863, 369, 115, 1, 62, "Subtitle"],
Cell[12981, 372, 249, 3, 35, "Subsubtitle"],
Cell[13233, 377, 330, 7, 53, "Input"],
Cell[CellGroupData[{
Cell[13588, 388, 199, 3, 35, "Subsubtitle"],
Cell[13790, 393, 359, 8, 53, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14186, 406, 115, 1, 35, "Subsubtitle"],
Cell[14304, 409, 444, 10, 58, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14785, 424, 252, 3, 35, "Subsubtitle"],
Cell[15040, 429, 630, 16, 53, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15707, 450, 315, 4, 35, "Subsubtitle"],
Cell[16025, 456, 1204, 21, 58, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17266, 482, 350, 5, 35, "Subsubtitle"],
Cell[17619, 489, 1010, 18, 35, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18666, 512, 371, 5, 35, "Subsubtitle"],
Cell[19040, 519, 1012, 18, 35, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[20101, 543, 155, 2, 62, "Subtitle"],
Cell[20259, 547, 240, 3, 35, "Subsubtitle"],
Cell[20502, 552, 401, 8, 53, "Input"],
Cell[CellGroupData[{
Cell[20928, 564, 199, 3, 35, "Subsubtitle"],
Cell[21130, 569, 448, 9, 53, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21615, 583, 104, 1, 35, "Subsubtitle"],
Cell[21722, 586, 347, 9, 58, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22106, 600, 293, 4, 35, "Subsubtitle"],
Cell[22402, 606, 764, 16, 53, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[23215, 628, 139, 1, 62, "Subtitle"],
Cell[23357, 631, 175, 2, 35, "Subsubtitle"],
Cell[23535, 635, 442, 15, 35, "Input"],
Cell[23980, 652, 1602, 47, 64, "Input"],
Cell[25585, 701, 1983, 47, 212, "Input"],
Cell[CellGroupData[{
Cell[27593, 752, 57, 0, 35, "Subsubtitle"],
Cell[27653, 754, 1384, 45, 88, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29074, 804, 97, 1, 35, "Subsubtitle"],
Cell[29174, 807, 1254, 33, 135, "Input"],
Cell[30431, 842, 1184, 32, 111, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[31652, 879, 115, 1, 35, "Subsubtitle"],
Cell[31770, 882, 2134, 54, 327, "Input"],
Cell[33907, 938, 2166, 55, 349, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[36122, 999, 89, 1, 62, "Subtitle"],
Cell[36214, 1002, 908, 20, 83, "Input"],
Cell[CellGroupData[{
Cell[37147, 1026, 280, 7, 53, "Input"],
Cell[37430, 1035, 163, 2, 35, "Output"]
}, Open  ]],
Cell[37608, 1040, 222, 7, 56, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[37867, 1052, 73, 1, 62, "Subtitle"],
Cell[37943, 1055, 39, 0, 35, "Subsubtitle"],
Cell[37985, 1057, 573, 11, 35, "Input"],
Cell[CellGroupData[{
Cell[38583, 1072, 102, 1, 35, "Subsubtitle"],
Cell[CellGroupData[{
Cell[38710, 1077, 828, 19, 58, "Input"],
Cell[39541, 1098, 151, 2, 35, "Output"],
Cell[39695, 1102, 153, 2, 35, "Output"]
}, Open  ]],
Cell[39863, 1107, 249, 6, 53, "Input"],
Cell[CellGroupData[{
Cell[40137, 1117, 169, 3, 35, "Input"],
Cell[40309, 1122, 162, 2, 35, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[40520, 1130, 90, 1, 35, "Subsubtitle"],
Cell[CellGroupData[{
Cell[40635, 1135, 814, 21, 99, "Input"],
Cell[41452, 1158, 149, 2, 35, "Output"],
Cell[41604, 1162, 149, 2, 35, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[41790, 1169, 1655, 46, 167, "Input"],
Cell[43448, 1217, 149, 2, 35, "Output"],
Cell[43600, 1221, 151, 2, 35, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[43788, 1228, 1480, 31, 35, "Input"],
Cell[45271, 1261, 411, 8, 58, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[45719, 1274, 1482, 30, 35, "Input"],
Cell[47204, 1306, 355, 7, 35, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[47608, 1319, 102, 1, 35, "Subsubtitle"],
Cell[47713, 1322, 418, 11, 35, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[48180, 1339, 112, 1, 62, "Subtitle"],
Cell[48295, 1342, 6012, 148, 399, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[54344, 1495, 185, 3, 62, "Subtitle"],
Cell[54532, 1500, 2031, 40, 190, "Input"],
Cell[CellGroupData[{
Cell[56588, 1544, 2404, 60, 164, "Input"],
Cell[58995, 1606, 278, 6, 27, "Print"]
}, Open  ]],
Cell[CellGroupData[{
Cell[59310, 1617, 1781, 47, 102, "Input"],
Cell[61094, 1666, 5338, 91, 302, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[66469, 1762, 218, 5, 35, "Input"],
Cell[66690, 1769, 133, 2, 35, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
